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Abstract 

A set of equations is developed to describe a curve in space given the 
curvature re and the angle of rotation 9 of the osculating plane. The set 
of equations has a solution (in terms of re and 9) that indirectly solves the 
Frenet-Serret equations, with a unique value of 9 for each specified value 
of r. Explicit solutions can be generated for constant 9. The equations 
break down when the tangent vector aligns to one of the unit coordinate 
vectors, requiring a reorientation of the local coordinate system. 

1 Introduction 

Given the curvature re and torsion t, the Frenet-Serret equations 1 describe a 
curve in space parameterized by the arc length s: 



dT 

ds 
dN 

ds 
dB 

ds 
dR 

ds 



kN; (1) 

-kT+tB; (2) 

-tAT; (3) 

T. (4) 



Here R, T, TV, and B are the position, tangent, normal, and binormal 
vectors, respectively. These equations have no explicit solution (in terms of re 
and r) for the general case, although solutions for special cases exist 2 . 

It is shown here that a set of equations can be developed to describe a curve 
in space given the curvature re and the angle of rotation 9 of the osculating 
plane. The set of equations has a solution (in terms of re and 9) that indirectly 
solves the Frenet-Serret equations, and has a unique 9 for every value of r. 
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Figure 1: Local coordinate system; T is normal to the plane containing j & k' . 



Many problems 3-7 involve the use of the Frenet-Serret equations, requiring 
numerical approximations or the use of helical arc segments (each having con- 
stant t and n). Specifying k and 9 to generate a solution may be useful if r is 
not initially known. The torsion r can then be determined from 9. 

2 Mathematical Development 

A local coordinate system having the property T = % (figure 1) supports the 
definition of N: 

N = j' cos 9 + k 1 sin 9. (5) 

The curvature k and the angle of rotation 9 of the osculating plane (con- 
taining N and T) characterize the curve. When the plane containing T and 
the global coordinate j is normal to k! (figure 2) then 
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Figure 2: Angular orientation of the local coordinate system with respect to the 
global coordinate system. 



k = irw, Z ,-| = / = ■ (6) 



Equation (6) breaks down when T — ±j, requiring an alternate expression 
for k' (developed in section 4). However, when T ^ ±j, 

Substituting (7) and (6) into (5): 

ldTi -T fc sin 6-TiTj cos 9 

N i - --jt / ; W 



K d-S _ rp2 



1 dT. 



Nj = -^=cose^-T- (9) 

ldT fc TismO-TjTkCose 
Nk = — i— = , • (10) 

k as A — ^ 2 
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Equation (9) can be integrated directly: 



leading to 



k cos 9 da: 



cos 9 da; 



sin 1 Tj = sin 1 T j0 + f k 

J Sf) 

J Sn 



Tj = sin 



cos 9da 



= sin 8: 



Tj = Tj cos / k cos 6da + A / 1 — T 2 n sin / k cos 6da. 



L j0 ' 



Equation (8) is solved by noting that T k = y 1 - Tf - Tf = ^cos 2 
and introducing the variable (3 so that 



Ti = cos 8 cos /?; 
Tfe = cos 5 sin p\ 



Substituting into (8): 



—— = — KcosCsind cos p — coso snip— 
as as 

— k sin cos 5 sin /3 — /c cos cos <5 cos /3 sin 8 
cos <5 

= — k sin sin f3 — k cos cos /? sin <5. 



Equation (15) simplifies to 



d(3 k sin 9 
ds cos 5 ' 



or 



= Po + / pCicr = COS , , 

J cost) \coso/ 



so that 



cos 5 f s Ks'm6 cos 8 /" Ksmt 
Ti = Ti n — cos / —da — T^g —sin / —da; 



' cos Sq 



cos 8 



COS(5n 



cos 8 



where 



COS 



5 = Jl — Tj cos / k cos 9da — Tj sin / k cos 9da. 

J sn J sn 



(19) 



The solution for Tk follows from (14) and (16): 
, cos 6 



«o s: 
COS Oi 



5 f s ksi 

- COS / 

'0 Js CO! 



sin 9 cos S 

da + XL — sin 



cos S 



' COS Sn 



r s Ksin6» 
J sa cos 8 



da. 



'so — —°0 

It can be easily verified that (12), (18), and (20) meet the requirement: 

dT 



(20) 



ds 



(21) 



Generating an expression for the torsion r requires first computing N by 
substituting (12)-(14) into (8)-(10): 



Ni 



cos 9 sin S cos (3 — sin (3 sin 0; 



Nj — cos cos 5; 

Nk — — cos sin 5 sin /3 + cos /3 sin 6. 
Next, B — T X N: 



(22) 
(23) 
(24) 



.Bj = sin S sin cos /3 — cos 9 sin /?; 

Bj = — cos (5 sin 9; 

Bk = sin 5 sin sin /3 + cos cos [3. 

Equation (28) expresses the torsion as a function of 9: 



T = 



dB 



ds 



d9 

= — k tan o sin ( 

ds 



Equation (29) expresses r in terms of components of T and B. 



d9 
ds 



1-Tf 



(25) 
(26) 
(27) 



(28) 



(29) 



Finally, (2) serves as a check on the solutions for T, TV, B, and r. 



3 Discussion 



Integrating (29) leads to the following expression for 9: 



9 = 9 



kT j b j 



da. 



(30) 
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Equation (30) indicates a unique value of 9 for each specified value of t when 
T 3 ^ ±1. Thus, (12), (18), and (20) indirectly solve (l)-(3). 

The angle 9 can also be expressed in terms of components of T, TV, B: 



sin 



-i 



— -=^= = cos — -=^= 



tan- (31) 



3.1 Constant 9 



An explicit solution often results when 9 is constant. Setting Tio = 1, so that 
00 = do = (and setting sq = 0) leads to 



S= k((j) cos 9oda; 
Jo 

(i = 2tan6>otanir 1 (tan (5/2) ; 



so that 



Tj = cos [2 tan 9 tanh 1 (tan (5/2)] cos / /c(er) cos 9 da; 

Jo 

Tj = sin / n(a) cos #od<r; 
Jo 

Tk = sin [2tan6*otanh _1 (tan 5/2)] cos / n(a) cos9oda. 

Jo 



The torsion becomes 



t(s) = — k(s) sin 6*o tan / k(ct) cos 9od<r. 
Jo 

As an example, when 



k = Koe , 



(32) 
(33) 

(34) 
(35) 
(36) 

(37) 
(38) 



Ti = cos 

Tj = sin 
Tk = cos 



2 
2 



erf(s) cos # 

erf(s) cos #o 
erf(s) cos 9q 



cos 



2 tan O tanh 1 ( tan 



erf(s) cos 9 



(39) 
(40) 



sin 



t(s) = — «o e s sin #o tan 



2 tan 6*o tanh 1 ( tan 



erf(s) cos #o 



erf(s) cos 9q 



(41) 
(42) 



G 



3.2 Constant k 



When k = kq but 9 ^ 9 , the solution will typically involve undetermined 
integrals. For example, when k = Kg and 9 — kqs, 



and 



Ti = cos (sin kqs 



) cos / 
Jo 



k sin n a 
cos(sinKocr) 



da; 



Tj = sin (sin kos) ; 

rp , ■ \ ■ f S k sin n Q a 

Tk = cos (sinKos) sin / — -da; 

J Q cos (sin Koa) 



t(s) = k — no tan (sin k s) sin k s. 
3.3 Constant k and 9 



(43) 
(44) 
(45) 

(46) 



When k = k q and 6 = 9 , (34)-(37) become: 

Ti = cos (koscos6>o) cos 2tan$o tanh -1 ^tan 
Tj = sin (k s cos 9 ) ; 

T k = cos (k scos9 ) sin 2tan# tanh~ 1 ^tan 



Hps 

KqS 



cos 9q 



cost 



(47) 
(48) 
(49) 



t(s) = — kq sin 9$ tan (kqs cos 9q) . (50) 

When #o = 7r/2, t(s) = 0, confining T and AT to a plane. When T aligns 
with j, t — > oo in (50), and the equations break down. 

4 Alternate Set of Equations 

The equations break down when Tj — > ±1, requiring a different orientation for 
the local coordinate system. The angle of rotation of the osculating plane is 
designated <p here. In general, 4> ^ 9, reflecting differences in angular orientation 
between the local and global coordinate systems for the two cases. Defining k! 
as the normal to the plane containing T and i, i.e., 



k' 



i x T -jT k + kTj 
\i x T| - v/I - 7f ' 



The j' unit vector becomes: 



k' xT = 



—i (l — Tf) + jTjTj + kTjT k 



(51) 



(52) 
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Substituting into the expression for N: 



at - l^El - -T k sin c^ + TiTj cos c/) 
j k ds Vl=T? ' 

ldT k Tj sin <j> + T{T k cos <f> 

n * = k-*t = — 7r=w — ■ (55) 

Equations (53)-(55) have the following solution: 

T 4 -sin 7 ; (56) 

Tj = cos 7 cos a; (57) 

Tfc = cos 7 sin a; (58) 

Ni = — cos 7 cos 0; (59) 

Nj — sin 7 cos a cos — sin a sin 0; (60) 

Nk — sin 7 sin a cos + cos a sin 0; (61) 

Bi = cos 7 sin 0; (62) 

Bj = sin 7 cos a sin <p + sin a cos <fi; (63) 

= — sin 7 sin a sin <j> + cos a cos (f>; (64) 



Here 



_- K tan7si n 0=-- T - w . (65) 



= sin 1 T i0 — J k cos 4>da; (66) 



cos" 1 ' 3 



's 

Tj \ f s k sin <fi 

= cuo + / «cr; (67) 

cos 7/ J So cos 7 

m m cos 7 /" s Ksin</> , „ cos 7 f s Ksmcb , , „ 

Tj = T jn -cos / -da-T kn -sin / -dcr; (68) 

cos 70 J so cos 7 cos 70 J so cos 7 

0^-tan- 1 ^. (69) 

Even though 9 and <f> both represent the angle of rotation of the osculating 
plane, (31) and (69) differ because of differences in angular orientation of the 
local coordinate system. 

When Tj — ► ±1 or Tj — > ±1, switching from one set of equations to another 
avoids numerical difficulties. 
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5 Concluding Remarks 



Unlike the Frenet-Serret equations, (8)-(10) are nonlinear, and do not involve 
N, _B, or r. The solution (in terms of n and 9) indirectly solves the Frenet- 
Serret equations, and leads to a precise definition of t as a function of k and 9. 
A unique value of 9 can be obtained for each specified value of r through a first 
order ordinary differential equation. The equations break down when T — > ±j, 
requiring an alternative set of equations that break down when T — > ±i. The 
expressions for the angle of the osculating plane in the two approaches differ 
because of differences in the angular orientation of the local coordinate system. 
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